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ABSTRACT 


Measuring redshifted CO line emission is an unambiguous method for obtaining an ac¬ 
curate redshift and total cold gas content of optically faint, dusty starburst systems. Here, we 
report the first successful spectroscopic redshift determination of AzTEC J095942.9H-022938 
(“COSMOS AzTEC-1”), the brightest 1.1mm continuum source found in the AzTEC/JCMT 
survey (Scott et al. |2008| l, through a clear detection of the redshifted CO (4-3) and CO (5- 
4) lines using the Redshift Search Receiver on the Large Millimeter Telescope. The CO 
redshift of z = 4.3420 ± 0.0004 is confirmed by the detection of the redshifted 158 fim 
[C II] line using the Submillimeter Array. The new redshift and Herschel photometry yield 
Lfir = (1-1 ± 0.1) X IO^^Lq and SFR « 1300Mq yr“^. Its molecular gas mass de¬ 
rived using the ULIRG conversion factor is 1.4 ± 0.2 x IO^^Mq while the total ISM mass 
derived from the 1.1mm dust continuum is 3.7 ± 0.7 x assuming Td = 35 K. Our 

dynamical mass analysis suggests that the compact gas disk (r « 1.1 kpc, inferred from dust 
continuum and SED analysis) has to be nearly face-on, providing a natural explanation for 
the uncommonly bright, compact stellar light seen by the HST. The [C II] line luminosity 
L[cii] = 7.8 ± 1.1 X 10® Lq is remarkably high, but it is only 0.04 per cent of the total 
IR luminosity. AzTEC COSMOS-1 and other high redshift sources with a spatially resolved 
size extend the tight trend seen between [C II]/EIR ratio and S fir among IR-bright galaxies 
reported by Dlaz-Santos et al. ( 2013[ ) by more than an order of magnitude, supporting the 
explanation that the higher intensity of the IR radiation field is responsible for the “[C 11] 
deficiency” seen among luminous starburst galaxies. 

Key words: galaxies: high-redshift - galaxies: starburst - galaxies: distances and redshifts - 
galaxies: individual: AzTEC J095942.9H-022938 - sub millimetre: galaxies - radio lines: ISM 


1 INTRODUCTION 
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Recent studies of cosmic star formation history and galaxy mass 
build-up have shown a remarkably tight correlation between star 
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formation rate (SFR) and stellar mass (M*), also known as star 
formation “main sequence”, for galaxies with M* up to 
extending out to 2 ~ 6 (see |Steinhardt et al.|2014[ [Salmon et al.| 
|2015| and references therein). A substantial population of quiescent 
galaxies with M, ^ are also found to 2 ; ~ 4, suggest¬ 

ing rapid formation and quenching of massive galaxies at 2 ~ 6 
or earlier ( [Whitaker et al.|20r3l jStraatman et al.|20T4| l. Given the 
constraints on rapid formation and cessation of stellar mass build¬ 
up and their compact morphology, intense starbursts and feedback 
driven by a rapid gas accretion are thought to be important in this 
process (see jWilliams et al.|2014| and references therein). 

The submillimetre galaxies (SMGs) are natural laboratories 
for testing this hypothesis and probing the details of the physical 
processes that govern this rapid build-up and quenching of massive 
galaxies. SMGs are identified by their large FIR luminosity, which 
is widely interpreted to b e powered by intense star formation with 
SFR > yr"^ jBlain et al.|2^|Yun et al.|2012| |Kirk- 


jpatrick et al.|2012T l. Wide area surveys by the Spitzer and Herschel 
Space Telescopes have shown that these luminous IR galaxies ac¬ 
count for a significant fraction of the Cosmic IR background (|Pen-| 
jner et al.|20I l[[Bethermin et al.|201^ , suggesting that they are an 
important component of the cosmic mass build-up history at 2 ^ 1 
( |Le Floc’h et~aL][2005] jCaputi et al.|[2007i [Magnelli et tlllTOTT] ). 

Because of their faintness in the optical bands, their precise red- 
shift distribution is poorly determined (see below), but [Toft et al.j 
( [2014^ have found that the well established population of massive 
(M, > lO^^M©) compact quiescent galaxies at 2 ~ 2 can be fully 
accounted by the known SMGs at 3 < 2 < 6 in terms of their 
abundance and stellar population, and jSimpson et al.] ( |2014^ repro¬ 
duce the local elliptical luminosity function by passively evolving 
the population of bright SMGs. 

The study of the brightest SMG found in the COSMOS (|Scov-| 
jille et al.||2507| ( field, AzTEC 1095942.9-1-022938 (“COSMOS 


AzTEC-1 ” hereafter), exemplifies the major challenge behind mak¬ 
ing this important connection between SMGs and the rapid build¬ 
up of stellar mass in galaxies. Eirst discovered by the AzTEC 
COSMOS survey using the James Clerk Maxwell Telescope (|Scott| 
[er^[2008) , COSMOS AzTEC-1 is one of the brightest SMGs 
known and is particularly well studied because of the extensive 
deep multi-wavelength data readily available in the COSMOS field 
(see jSmolcic et al.|[2011| l. Unlike many other SMGs (including 
those discovered by Herschel) that suffer from low angular reso¬ 
lution of single dish telescopes and source blending, the location 
of this SMG is known to better than 0.1” accuracy because of a 
dedicated interferometric imaging survey done using the Submil¬ 
limeter Array (SMA) by jYounger etar] ( |2007l l. COSMOS AzTEC- 

I is the only object among the 7 AzTEC sources imaged with SMA 
by Younger et al. that has an unambiguous optical counterpart, and 
this relatively bright (m — 25.3 mag [AB] in the HST F814W 
band) source is extremely compact, ~0.2" (~1.5 kpc) in diame¬ 
ter, comparable to the sizes of the massive compact galaxies found 
at 2 = 2 ~ 4. A 4 hour long exposure with DEIMOS on Keck 

II telescope by jSmolcic et al.j j201 1| | did not yield any emission 
lines, and the continuum break near 6700 A was interpreted as the 
blue cutoff of Ly-a at 2 = 4.65. Also, using 31 NUV-NIR photo¬ 
metric measurements, Smolcic et al. derived a photometric redshift 
of 2 = 4.64 with a secondary peak at 2 = 4.44. Their attempt 
to confirm this redshift by CO spectroscopy using CARMA and 
PdBI interferometers failed to detect any CO emission in the red¬ 
shift range of 4.56 < 2 < 4.76 and 4.94 < 2 < 5.02. Later, jlon^ 
jet al.| ( [20T^ expanded the CO line search using the Nobeyama 45- 
m Telescope to 4.38 < 2 < 4.56, but they also failed to detect 


a CO line. Therefore, the redshift of this arguably the best studied 
AzTEC SMG in the COSMOS field still lacks a spectroscopic con¬ 
firmation despite nearly 10 years of efforts using some of the most 
powerful astronomical facilities available. With the exception of the 
SMG COSMOS AzTEC-3, which is recently shown to be part of 
a large scale structure at 2 = 5.3 jRiechers et al.|201^|2014^ , the 
situation is essentially the same for the remaining AzTEC sources 
imaged by the SMA as well - all are expected to be at 2 ^ 3 be¬ 
cause of their faintness in the optical and the radio bands, with a 
much worse prospect of yielding a spectroscopic redshift. 

In this paper, we report the first successful spectroscopic red¬ 
shift determination of COSMOS AzTEC-1 obtained using the Red¬ 
shift Search Receiver (RSR) on the Large Millimeter Telescope Al¬ 
fonso Serrano (LMT), which is a ultra-wide bandwidth spectrome¬ 
ter designed to conduct a blind search for redshifted CO lines from 
molecular gas rich galaxies. We also report the confirmation of the 
CO redshift through the detection of the redshifted 158 pm [C II] 
line using the Submillimeter Array. We interpret the observed CO 
and [C II] luminosity in terms of a highly concentrated and in¬ 
tense starburst, fueled by the CO and [C II] emitting gas and its 
properties, including the “[C II] deficiency” in COSMOS AzTEC- 
1. Throughout this paper, we assume flat ACDM cosmology with 
Om = 0.3 and J Tq = 70 km s~^ Mpc“^ and Kroupa initial mass 


function (IME; Kroupa 


2001 
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2 OBSERVATIONS 

2.1 Redshift Search Receiver Observations on the LMT 

The Redshift Search Receiver observations of COSMOS AzTEC-1 
were conducted in January and February 2014 as part of the Early 
Science program at the Large Millimeter Telescope ( [Hughes et al.j 
[2010[ l. The Redshift Search Receiver (RSR) consists of two dual 
polarization front end receivers that are chopped between the ON 
and OFF source positions separated by 76” in Azimuth at 1 kHz 
rate using a ferrite switch, producing a flat baseline over the entire 
38 GHz (73-111 GHz) bandwidth and always integrating on-source 
(see [Erickson et al.|20(37{ [Chung et al.[[2009| for further descrip¬ 
tions of the instrument). The ultra-wideband backend spectrometer 
covers the entire frequency range between 73 and 111 GHz simul¬ 
taneously with 31 MHz {R — 3000 or 100 km s“^ at 93 GHz) 
spectral resolution. During this Early Science phase operation, only 
the inner 32 metre diameter section of the telescope surface is illu¬ 
minated, leading to an effective beam size of 20” at 110 GHz and 
28” at 75 GHz. A total of 290 minutes of on-source integrations 
were obtained over 3 different nights, mostly in excellent weather 
with Tsya « 90 K (t 225 GHz = 0.05 — 0.1). Telescope pointing 
was checked every 60-90 minutes by observing the nearby QSO 
J0909-H013. 

Data were reduced and calibrated using DREAMPY (Data 
REduction and Analysis Methods in PYthon), which is the RSR 
data reduction pipeline software written by G. Narayanan. After 
flagging any data adversely affected by a hardware or software 
problem, a linear baseline is removed from each spectrum. The 
final spectrum shown in Figure [T] was obtained by averaging all 
spectra using the l/a^ weight, and the resulting final rms noise is 
a = 0.13 mK. The measured gain as a function of elevation and 
frequency using Uranus and MWC349A is flat, 7 Jy/K (in TX unit) 


^ Stellar mass in Kroupa IMF is 38 per cent smaller than that of the Salpeter 
IMF-i.e., Mt,{Kroupa) = 0.62M*(Salpeter). 
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between the elevation range of 30-75 degree, where all observations 
were made. 


2.2 Submillimeter Array 

Spectroscopic imaging observations of COSMOS AzTEC-1 were 
obtained on 1 March 2014 using the using the Submillimeter Array 
(SMA; |Ho et al.|2004| >, the 8-element interferometer located near 
the summit of Mauna Kea, Hawaii. The SMA was in a close pack 
configuration with projected baselines ranging from 6 to 45 metre 
(mean ~21 metre). The array was operated using two orthogonally 
polarized SIS receivers on each antenna, each tuned to 355.7 GHz 
within the 2 GHz wide upper sideband, the expected frequency for 
the redshifted [C II] line based on the results of the LMT observa¬ 
tions. The lower sideband was also captured allowing a sensitive 
measure of the thermal continuum near 345.4 GHz (1850 GHz = 
162 /rm in the source frame). The raw spectral resolution was 3.25 
MHz uniform over both sidebands, or about 2.75 km s“^ around 
the [C II] line. 

The spectral response was calibrated using observations of the 
bright QSO 3C 84, and the flux density scale was calibrated from 
measurements of the Jovian moon Callisto, known to within 5 per 
cent in the submillimetre bands (based upon SMA observations, see 
ALMA Memo 59^. Observations of the target were interleaved 
with measurements of QSOs J0909-I-013 (0.326 Jy) and J1058-I-015 
(2.31 Jy) for use in calibrating the complex gains due to instrumen¬ 
tal and atmospheric effects. The observations were obtained in very 
good weather, with T 225 GHz ~ 0.075. The total on-source integra¬ 
tion time was 5.2 hours. The complex visibility data were calibrated 
within the MIR reduction packa g^ and the calibr ated visibilities 
were then exported to MIRIAD ( |Sault et al.|1995| l for resampling 
to a common spectral grid. 

The continuum and the [C II] spectral line image cube 
were produced using the Astronomical Image Processing System 
(AlPSj^task IMAGE. Since no spatial details are expected to be 
revealed at the angular resolutions achieved (see below), natural 
weighting was used in the mapping in order to maximize sensitiv¬ 
ity. The 345 GHz continuum image with an effective bandwidth of 
2 GHz has a synthesized beam of 5.8” x 3.4" (PA = 60°) and Icr 
noise of 1.4 mjy beam“^. The spectral resolution of the line data 
is 11.865 MHz (10 km s"^ at 355 GHz). The final [C II] spectral 
line cube was produced by first subtracting the continuum from the 
uv-data and then averaging over 8 channels with an increment of 
4 channels, covering the frequency range between 354.4347 GHz 
and 356.2381 GHz. The resulting cube has a spectral resolution of 
94.92 MHz with a synthesized beam of 5.7" x 3.3" {PA = 59°) 
and la noise of 4.2 mJy beam“^ in each channel. 


3 RESULTS 

3.1 CO Redshift and Line Luminosity 

The final RSR spectrum of COSMOS AzTEC-1 shown in Figure[T] 
has two emission lines clearly above the noise level. A straightfor¬ 
ward interpretation of the spectrum is that these are two redshifted, 
adjacent rotational transitions of CO, and the separation between 


^ http;//library.nrao.edu/public/memos/alma/ 
memo594.pdf 

https : / / WWW . cf a . harvard. edu / -cqi/mircook. html 
^ http://WWW.alps.nrao.edu/index.shtml 
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Figure 1. RSR spectrum of COSMOS AzTEC-1. The two spectral features 
well above the noise level are interpreted as CO (4—3) and (5^) lines at 
2 = 4.342. The redshifted 492 GHz [C I] line should appear at 92.13 GHz 
(marked with an arrow). 
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Figure 2. Zoom-in views of the two CO lines in the RSR spectrum of COS¬ 
MOS AzTEC-1. The x-axis is in velocity offset with respect to the system¬ 
atic redshift of 2 = 4.3420 (vertical long-dashed lines). 


the two lines, Ap = 21.565 GHz, corresponds to the expected fre¬ 
quency offset between CO J = 4 —> 3 and J = 5 —>■ 4 transitions 
at 2 = 4.342 (see Appendix for a detailed discussion on the red¬ 
shift determination from an RSR spectrum). Both lines are fully re¬ 
solved by the RSR (see Figure 2 1 , and the best fit Gaussian param¬ 
eters are summarized in Table 1 The CO (4-3) line is centered at 
V = 86.3085 GHz, corresponding to a redshift of 4.3418 ± 0.0006 
while the CO (5^) line is centered at tz = 107.8739 GHz, at a 
redshift of 4.3421 ± 0.0006. The best-fit linewidths (FWHM) are 
380 km s“^ and 364 km s“^, respectively, in good agreement as 
expected if they are two CO transitions from the same galaxy. The 
redshifted 492 GHz [C I] line, which should appear at 92.130 GHz, 
is undetected with S'jc i\/Sco{ 4 ,- 3 ) ^ 0.45, in line with the mea- 
I] line strengths in other high redshift galaxies jWalter et| 
•^[C I]/'S'cO(3-2) ~ 0.3). 

an ultra-wide spectrum produced by the RSR, the redshift 
information is also present in weak lines such as [C I] that are not 
formally detected individually as well as in bright lines such as 
CO. We have developed a method to exploit all spectral information 
present in the RSR data by cross-correlating the observed spectrum 
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Table 1. Summary of the CO and [C II] Line Measurements 


Line 

"CO/[C II] 
(GHz) 

^CO/[C II\ 

AC 

(km s“®) 

5Ay 

(Jy km s~^) 

Leo /[Oil] 
(10®L©) 

r / 

^CO/[C II] 

(10^^ Kkms-1 pc^) 

^H2 

(10®°M©) 

CO (4-3) 

86.3085 

4.3418 ± 0.0006 

380 

1.75“ ± 0.24 

2.5 ±0.3 

7.8 ± 1.1 

14 ±2 

CO (5^) 

107.8739 

4.3421 ± 0.0006 

364 

1.55“ ± 0.22 

2.8 ± 0.4 

4.4 ± 0.6 

8.8 ± 1.2 

[CII] 

355.8038 

4.3415 ± 0.0003 

366 

13.05 ±0.70 

78 ±5 

3.6 ±0.2 

- 


(a) A conversion of 7 Jy/K is adopted to convert the measured antenna temperature in Ta* to flux density, using the calibration factor derived between 

December 2013 and January 2014. 

, and Mh 2 is derived using the “ULIRG” conversion factor 


(b) is estimated using the average line ratios for SMGs Carilli & Walter 


2013 


o-co = O.SMq (K km s ^ pc^) ^ — see§[4G 


with a theoretical or an empirical spectral template (see |Yun et al.| 
|2007| l. A detailed analysis of the cross-correlation amplitude along 
with the expected CO line multiplicity and the redshift constraints 
from the radio-millimetric photometric redshift analysis uniquely 
identifies the 2 = 4.342 solution with a total S/N — 9.0 (see 
Appendix). As shown in the zoom-in inset of Figure [AT| this red¬ 
shift peak is well-resolved by the cross-correlation analysis with 
a spread in redshift between 4.338 and 4.347 (FWHM). It is well 
centered on the redshift of AzTEC-1 derived from fitting the in¬ 
dividual CO lines, 2 = 4.3420 ± 0.0004 (see Table [TJ, but the 
width of the distribution is nearly 10 times larger than the uncer¬ 
tainty from the individual line fitting, indicating that the width of 
the cross-correlation amplitude arises from the finite width of the 
CO lines (~375 km s“^) rather than reflecting the uncertainty in 
the redshift determination. 

The CO line redshift of 2 = 4.3420 ± 0.0004 we derive 
from the RSR spectrum is significantly lower than the Lyman-a 
break based redshift of 4.650 ± 0.005 o r the optical/IR photo met- 
ric redshift of 2 = 4.64('lQ Qg reported by Smolcic et al. (2011 1 , and 


naturally explains why their CARMA and PdBI CO line searches 
failed. Their photometric redshift analysis produced a secondary 
solution at 2 = 4.44, which is much closer to our CO redshift. 
Although |Iono et ah] ^2012) had the right idea to search for a 
CO line near this secondary redshift peak, they missed detecting 
the CO (5-4) line just outside their search range. In either case, 
the redshift adopted by |Smolcic et al.| )201 1) is close enough to 
the actual CO redshift that their analysis of stellar mass and IR 
spectral energy distribution is still mostly valid, and their conclu¬ 
sion that AzTEC-1 is an extremely young (^ 50 Myr), massive 
(Mt ~ IO^'^Mq), and compact (^ 2 kpc) galaxy with a star for¬ 
mation rate of SFR ~ 10®Mq yr“® still holds. 

The CO line luminosity Leo in T© can be computed from 
the measured line integrals and the CO redshift using Eq. (1) by 
[Solomon et al.| ( [i997) as, 

Leo = 1-04 X 10 ^ Sco^Vi'o{l + z) ^ D\ [L©] (1) 

where Sco^V is the measured CO line integral in Jy km s“®, uq 
is the rest frequency of the CO transition in GHz, and Dl is the 
luminosity distance in Mpc. As summarized in Table [T] Leo is 
(2.5 ± 0.3) X 10®L© and (2.8 ± 0.4) x 1 O®L 0 for the CO (4-3) 
and CO (5-4) transitions, respectively. 

Total molecular gas mass is related to the quantity L'eo (see 
Eq. (3) by jSolomon et al.|1997) , 

L'eo = 3.25 x 10^5coAI/iz-(.® (1 -b 2 )“®!)! [TC fern s-®pc®]( 2 ) 


where fobs = + z) is, the observed line frequency in GHz. 

The derived CO luminosities are L 0 q( 4 _ 3 ) = (7.8 ± 1.1) x 10®® 
Kkms“® pc® andL 0 Q(.g_ 4 ) = (4.4±0.6) x 10®® Kkms“® pc®. 



Figure 3. SMA 345 GHz continuum image of COSMOS AzTEC-1 at 
5.8” x 3.4" resolution. Contours correspond to 2a, 4cr, her, 8(t, & lOtr 
(a = 1.4 mJy/beam). The HST i-band image is shown in greyscale. 


These quantities can be converted to molecular gas mass Mh 2 with 
several assumptions. Given the large uncertainties involved in this 
conversion, the total gas mass estimation from L'eo ** deferred to 
a discussion later (see § [4.1 [below). 


3.2 [C II] Line and 345 GHz Continuum 

While the two CO lines detected with the RSR on the LMT look 
quite good, confirming the CO redshift of 2 = 4.342 and rul¬ 
ing out the possibility of a chance superposition of two unrelated 
low redshift (2 < 3) CO sources along the same line of sight re¬ 
quires a confirmation with another emission line. The redshifted 
[C II] line (i/Q = 1900.537 GHz) falls in the middle of the SMA 
350 GHz receiver band, which also overlaps with the frequency 
coverage of the SMA 400 GHz receiver. The SMA spectrum ob¬ 
tained (Figure]^ shows a clearly detected and fully resolved [C II] 
line near 355.7 GHz with S/N ~ 15, confirming the CO red¬ 
shift. Three other redshifted [C II] line sources (the 2 = 4.7 QSO 


BR1202-0725 Pono et al.|2006) , the 2 = 5.24 lensed Herschel 


source HLSJ091828-1-514223 ([Rawle et al.[2014 

, and the 2 = 4.68 

lensed SMG HLS1-MACSJ2043 IZavala et al. 

|2015|) have been 


detected by the SMA before, and this new [C II] detection of COS- 
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Figure 4. A [C II] spectrum of COSMOS AzTEC-1 obtained using the 
SMA. The [C II] line redshift of z^c II] = 4.3415 ± 0.0003 and AV" = 
366 km s“^ are in excellent agreement with the RSR CO line measure¬ 
ments (see §|34j. 


MOS AzTEC-1 nicely demonstrates the intrinsic brightness of the 
[C II] line and the excellent sensitivity of the SMA for studying 
high redshift [C II] sources. 

Because the [C II] line is detected with a significantly higher 
S/N and a higher spectral resolution than the RSR CO observa¬ 
tions, this SMA [C II] spectrum can reveal much more than sim¬ 
ply confirming the redshift of COSMOS AzTEC-1. The [C II] line 
shown in Figure|^is slightly asymmetric, spanning 520 MHz (438 
km s“^). The center of the line at full-width-zero-intensity (FWZI) 
corresponds io v = 355.74 GHz oi z — 4.3425, which agrees 
very well with the CO redshift. The best fit Gaussian model for 
the line yields a mean redshift of Z[c ii] = 4.3415 ± 0.0003 with 
AH = 366 km s“^, reflecting the asymmetry with brighter emis¬ 
sion on the blue-shifted side of the line. The origin of this line 
asymmetry is unknown at the moment, and future high spectral res¬ 
olution measurements with a better S/N should reveal whether this 
asymmetry is also present in CO lines and potentially offer a useful 
insight into the spatial distribution between the CO and the [C II] 
emitting gas. 

The measured [C II] line integral of S[cii]^V = 13.05 ± 
0.70 Jy km s”'^ translates to [C II] line luminosity of L^c ii] = 
1.04 X IQ-^ Sic ii]AV uo{l + z)-^Dl = 7.8 ± 1.1 x IO^Lq. 
This is nearly 30 times larger than the line luminosity of both CO 
transitions detected with the RSR (see Table[^. To put this in per¬ 
spective, the [C II] line luminosity of COSMOS AzTEC-I is only 
a factor of 3 smaller than the total IR luminosity of an L* galaxy 
in the local universe ^Soifer et al.|1987HSaunders et al.|1990[|Yun| 
|et al.|2001^ . Indeed COSMOS AzTEC-I is extremely luminous in 
[C II] line, and this explains why the [CII] line is detected so easily 
by the SMA. Nevertheless, this [C II] line luminosity is only 0.04 
per cent of the total IR and bolometric luminosity (see below), and 
this “[C II] deficiency” is discussed further in detail in § |4.2| 

Both the [C II] line and the 345 GHz continuum detected are 
centered precisely on the position of the Hubble ACS i-band source 
as shown in Figure]^ The total measured 345 GHz continuum flux 
of 17.8 ±1.4 mjy is in excellent agreement with the 340 GHz con¬ 
tinuum fluxes of 15.6 ± 1.1 mJy previously reported by | Young^ 
|et al.| ( (2()07| >. This continuum is only marginally resolved by the 
longest baselines of the SMA with an inferred Gaussian source di¬ 


ameter of only 0.3”i Younger et al. 2008| . Both the line and contin¬ 
uum emission is unresolved by the 5.8" x 3.4" resolution of the 
new SMA data, as expected. 


3.3 IR Luminosity and SFR 

The extensive and deep multi-wavelength photometric data readily 
available in the COSMOS field allowed fSmolcic et al.| ( |20TTt to 
assemble an impressive array of spectral energy distribution (SED) 
data for their photometric redshift analysis as well as stellar pop¬ 
ulation modeling and bolometric luminosity estimation - see their 
Table 1 and Figure 4. By adding the Herschel SPIRE 250, 350, and 
500 fim photometry ( [Smith et al.|2012) as well as our new 345 GHz 
continuum measurement (see § |4.2^ , we have fully mapped the in¬ 
frared peak of the SED as shown in Figure]^ and a more reliable 
analysis of the dust heating and infrared luminosity is now possible. 

We model and interpret the observed SED using three differ¬ 
ent commonly used tools: a modified black body model, a starburst 
SED model byj Efstathiou et al.] (|2000[, and the GRASIL population 
synthesis and radiative transfer model ( [Silva et al.|1998| (. The modi¬ 
fied black body model characterizes only the far-IR part of the SED 
as dust processed radiation at equilibrium temperature weighted by 
luminosity. The two latter models aim to gain further insight into 
the nature of the luminosity sources by adding assumptions on the 
source geometry and star formation history. All these models are 
highly idealized, however, and these interpretations should be taken 
in the context of the assumptions adopted. 


3.3.1 Modified Black Body Model 

A common illustrative model for thermal dust emission from astro¬ 
nomical sources is modified black body radiation or “grey body” 
radiation. Following the classical derivation by [Hildebrand[ (198^ 
and adopting an emissivity function of the form Q{v) = 1 — 
(so that the emerging spectrum is pure black body 
spectrum aiv and Sv oc gj lower u) leads to a simple 

functional form of Td)[ 1 — e 2 ;p[—(^)'^], where 

Q.d is the solid angle of the source and B{v,Td) is Planck Function 
at frequency o and dust temperature Td ( [Yun et al.[2002[ . 

The best fit model that describes the observed SED between 
250 and 1000 fim photometry measurements is shown in Figure|^ 
characterized by dust temperature Td = 54 ± 3 K and emis¬ 
sivity /3 = 1.6 ± 0.2. The derived IR luminosities are Lm = 
1.4 X 10^® 1/0 and Lfir — 1.0 x IO’^^Lq, where Lm and Lfir 
are luminosities between A = 8 — 1000 fim and A = 40 — 120 fim, 
respectively. These luminosity estimates agree well with other esti¬ 
mates discussed below, although they are slightly smaller because 
this single dust component characterization does not account for 
the warm dust contribution in the mid-IR region. The derived dust 
emissivity/3 = 1.6 is similar to the commonly adopted value of 1.5 
and is known to be somewhat degenerate with Td in this formula¬ 
tion. The dust temperature of AzTEC-1 follows the general trend of 
increasing Td with Ljr reported by recent statistical studies such 
as by [Symeonidis et al.[ ( [2013') l and [Magnelli et ar]p012[[2014[ (. 
The derived dust temperature of 54 K for AzTEC-1 is higher than 
the average (Td) « 40 K for Lir = IO^^Lq sources at 2 = 2 
analyzed in these studies, and their predictions on the redshift evo¬ 
lution of Td with Lm differ slightly - sample selection is likely 
important. Far-IR data for 2 > 4 sources are rare because of the 
limitations of existing facilities, but our derived dust properties are 
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Figure 5. A spectral energy distribution (SED) of COSMOS AzTEC-1 from 
UV to radio wavelengths. Most of the photometry points shown are already 
summarized in Table 1 by |Smolcic et al.H2011) . The new 345 GHz SMA 
photometry and the Herschel SPIRE 250, 350, and 500 /im points (shown 
in red) are from the published catalog by |Smith et al.|j2012| clearly map 
out the dust peak, allowing us an accurate IR luminosity for the first time. 
The best fit modified black body model fitting the far-IR part of the SED 
is shown in magenta line (see the text for details). The “starburst” model 
SEDs by lEfstathiou et al.]j200Q^ with ages of 26, 37, 45, & 64 Myr are 
shown for comparison. 


similar to those of the seven z > 4 SMGs with L ~ ana- 

lyzed by lHuang et aT] ( [M^ Td = 40-80 K). 


3.3.2 Starburst SED Model by Efstathiou et al. (2000) 

Model SEDs of young stellar clusters embedded in a giant molecu¬ 
lar cloud by |Efstathiou et al^j2000[ > are shown in Figure]^ primar¬ 
ily for illustrative purposes and to compute the IR luminosity of 
COSMOS AzTEC-1. Although based on a relatively simple geom¬ 
etry and a highly idealized star formation history (a r-model), these 
“starburst” models as well as “cirrus” models with a lower opacity 
are shown to be remarkably effective in reproducing the observed 
SEDs of high redshift ultraluminous infrared galaxies (ULIRGs) 
and SMGs (see [Efstathiou & Siebenmorgen|200^ . The main im¬ 
pact of increasing starburst age is the build-up of the photospheric 
emission at wavelengths shorter than 3 /rm (with a corresponding 
increase in luminosity in the rest frame optical bands) and a system¬ 
atic shift of the dust peak to a longer wavelength as average opacity 
decreases and lower mass stars contribute more to the luminosity 
(i.e., cooler dust temperature). The success of these relatively sim¬ 
ple SED models can be attributed at least in part to the basic fact 
that the youngest stars dominate the luminosity and the detailed star 
formation history is largely washed out. Thus, we expect the IR lu¬ 
minosity and the current star formation rate to be reliable but the 
mass of the stars produced by the ongoing starburst to be less cer¬ 
tain. The SED with 45 Myr old starburst is in very good agreement 
with nearly every photometry data in Eigure|^ but this time scaling 
may be meaningful only in the context of this specific model. The 
radio to millimetre wavelength part of the SED is constructed using 
the well established radio-IR correlation for star forming galaxies 


as described by [Yun et al.H2002) , and the observed 1.4 GHz ra¬ 
dio emission is entirely consistent with the radio and IR luminosity 
being powered by a pure starburst. 

Using the best fit SED (“45 Myr”) model shown in Fig- 
ure[^ we estimate the luminosities of Lm = 1.5 x IO^^Lq and 
Lfir = 9.1 X lO^^I/Q, where Lm and Lf/a are luminosities be¬ 
tween A = 8 — 1000 /rm and A = 40 — 120 /^m, respectively. The 
star fo rmation rate deriv ed from Lm using the empirical calibra¬ 
tion by Kennicutt 1 1998 i.e., SFR = L/_r/(9.4x lO^L©) MQ/yr 


adjusted for Kroupa IMF) is 1596 MqIjt. Assuming z = 4.64, 
Smolcic et al. 1 2011 1 estimated Lm = 2.9 x lO'^^L©, nearly a fac¬ 


tor of 2 larger than our estimate when corrected for the new redshift, 
because they did not have the Herschel photometry to constrain the 
FIR peak. In comparison, the current SFR computed by the best 
fit SED model shown in Figure|^is 880 MqIjt, about a factor of 2 
smaller than estimated from the IR luminosity. The smaller current 
SFR of this model stems from the exponentially decreasing star 
formation scenario adopted by the model and may not be accurate. 


3.3.3 GRASIL SED Models 

GRASIL is a population synthesis code which predicts the SED 
of galaxies from far-UV to radio wavelengths jSilva et al.|1998| >. 
By allowing a wide range of geometry for gas/dust and stars and 
a realistic treatment of dust processing as well as a variety of 
star formation histories and IMFs, GRASIL can model SEDs of 
a wide range of plausible astrophysical scenarios. To model the 
observed SED of COSMOS AzTEC-1, we adopt a Schmidt type 
law (SFR{t) = uschMg^^) with an efficiency of uscH = 0.5 
Gyr~^ and an exponent k = 1 for the quiescent star formation 
history. The starburst component is modeled as an exponentially 
decreasing SFR{t) with an e-folding time tb, observed at different 
times (“ageb’’] after the outset of the burst. Both stellar and gas/dust 
sources are modeled with a King profile with core radii of r* and 
Cgas. Two different components of gas and dust are also consid¬ 
ered: (a) molecular clouds (MCs) where young stars are forming; 
and (b) “diffuse” or “cirrus” component that surrounds the MCs, 
old free stars, and exposed new stars with a large filling factor. A 
self-consistent radiative transfer calculation is performed to com¬ 
pute the emerging SED. Uncertainties in the model parameters are 
estimated from the analyses of 250 Monte Carlo realisations of the 
input photometry data. 

The best fit GRASIL starburst SED model for COSMOS 
AzTEC-1 is shown in Figure along with the photometry data 
used to constrain the model. An acceptable fit could be obtained 
without any AGN contribution, and the observed SED is consis¬ 
tent with the bulk of the luminosity being powered by a strong 
starburst, similar to the majority of local ultra-luminous infrared 
galaxies (ULIRGs; [Vega et al.|2008^ . The starburst component of 
the best fit model is characterized by an exponentially fading burst 
with an e-folding time of tb = 35 Myr, observed at agCb = 51 Myr. 
At these time scales, the majority of young stars are still embedded 
within their parent molecular clouds, and most of their luminosity 
emerges in the IR. Regardless of any detailed assumptions in the 
model, the global quantities such as luminosity and SFR should be 
quite robust. The total IR and FIR luminosities derived are Ljr = 
1.6±0.3 X mid Lfir = 1.2±0.2 x IO^^Lq, respectively, 

with the current star formation rate of SFR — 1320 ± 230 Mq 
yr“^. These IR luminosities are in good agreement with those de¬ 
rived in the previous section, and the model SFR is close to SFR 
derived from the IR luminosity, 1702 ± 296 Mq yr“^. 

There is little evidence to support any AGN activity in 
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Figure 6. Best fit GRASIL starburst SED model for COSMOS AzTEC-1 
is shown with a solid line. In addition to the photometry points shown in 
Fi gurej^ addition al photometry data in the UV and optical bands compiled 
by |Smolcic et al.|(2Ql 1) are also shown [in green]. The observed dust repro¬ 
cessed (FIR) light originates nearly equally from the diffuse medium (blue 
dashed line) and the molecular clouds (red dotted-dashed line). 


AzTEC-1, despite its large luminosity {Lm > While 

the light distribution seen in the HST i-band image (Fig.[^ is com¬ 
pact, it is clearly resolved with a diameter of 0.3^^ similar to other 
2 > 3 SMGs that are compact and clumpy with rg ^ 2 kpc l |Toft et| 
|al.|2014^ . The photometry data is sparse in the mid-IR (5-50 fim) 
range, and the presence of a heavily obscured AGN cannot he com¬ 
pletely ruled out hy this modeling. The red continuum between 0.8- 
8.0 fj.m can be interpreted as an indication of a buried AGN jLacy| 
let al.|2004| [Stern et al.|2005) , but it is also the characteristics of a 
heavily obscured young star clusters, commonly seen among most 
SMGs ( |Yun et al.|2008[ >. The observed radio continuum flux is en¬ 
tirely consistent with the expected supemovae rate (see Fig. 1 ^, and 
there is no room for any significant AGN contribution in the radio 
either. 

A good spectral coverage from the UV to the radio imposes 
strong constraints, not only on the global properties of the galax¬ 
ies but also on other important physical parameters. In the absence 
of an AGN, the fit to the near-IR and radio luminosities provides 
a strong constraint on the star formation rate, while the detailed 
shape of the SED is affected mainly by the value of the model pa¬ 
rameters such as star formation history and extinction. For instance, 
the SED shape in the UV range is determined by the geometry be¬ 
tween stars and dust while the SED in the optical range can help 
us to put constraints on the age of the old/intermediate age stellar 
populations. The optical depth mainly affects the mid-IR spectral 
range by varying the MC contribution. The best fit model core radii 
of the stellar and gas distributions are r* = 0.10 ± 0.02 kpc and 
rgas = 0.95 ± 0.42 kpc, respectively, in agreement with their sizes 
measured by the HST and SMA. The larger extent of gas and dust 
over the stars ensures an efficient obscuration of the stellar light, 
and the high mean opacity {Ay > 200 for the MC component), 
shapes the observed very red continuum SED in the rest frame op¬ 
tical and near-IR bands, as commonly seen in other high redshift 
SMGs (e.g., |Yun et al.|2008|[M^ . 

The total stellar mass derived by the GRASIL model M, = 
4.4 ± 0.7 X IO^^Mq is nearly twice as large as the estimate by 
ISmolcic et al.| ( |201 l| l, but this estimate depends strongly on the 
chosen star formation history and thus is not very secure. The total 
gas mass inferred from the GRASIL model is 3.6 ±0.6 x 10^^M©, 


and the 39 per cent of this total (1.4 x lO^^M©) is in the “dense” (or 
MC) phase directly fueling the star formation. As discussed below, 
deriving the total molecular gas mass from the new CO measure¬ 
ments requires several highly uncertain assumptions, and the gas 
mass estimate by the GRASIL model is near the high end of the es¬ 
timates derived using different methods and is close to the gas mass 
estimated from the d ust luminosity using the relation derived by 

of 


Scoville et al. 


1 2014 1 . The gas mass fraction fgas = 


45 per cent rs' sigmhcantly higher than the value derived in nearby 
galaxies of 10-20 per cent and is similar to the mass fraction found 
for the jz ~ 2 submillimetre galaxies (e.g., jTacconi et al.||2010| 
|Geach et al.|201 1) . 

About 50 per cent of the IR luminosity arises from the MC 
component while the other 50 per cent comes from the reprocessed 
light from free stars in the cirrus component. The average density 
of the dense MC component, umc = 7 x 10® cm“®, is higher 
than the critical density of CO (4-3) and CO (5-4) transitions, and 
this gas is capable of producing fully thermalized emission in these 
transitions. The GRASIL model does not constrain the nature of 
the cirrus component well, but the average density is expected to 
be 2-3 orders of magnitude lower when scaled by mass and sizes, 
and the two high J CO transitions are expected to be sub-thermally 
excited in this diffuse component. 

The best fit GRASIL mode for COSMOS AzTEC-1 is remark¬ 
ably similar to the best fit starburst SED models by |Efstathiou et al!] 
( |2000^ , as the gross shapes of the SED are driven largely by a young 
starburst embedded in dense gas clouds in both models. The best fit 
GRASIL model also requires a starburst history remarkably similar 
to the Efstathiou models discussed above, and this offers a plausi¬ 
ble explanation for the apparent success of the Efstathiou models 
despite its simplicity. With a more sophisticated parameterization, 
the GRASIL model can offer a more nuanced physical insight into 
the gas, dust, and stellar properties. 


4 DISCUSSION 

4.1 Gas Mass and Nature of the Molecular Gas Fueling the 
Luminosity 

A standard practice for deriving total molecular gas mass from 
a redshifted CO line measurement is to adopt a standard line 
ratio between different rotational transitions of CO to estimate 
the luminosity of the CO (1-0) transition Lco(r-o) *^hen 
to translate this line luminosity to molecular gas mass assuming 
an “aco” conversion factor (see a review by jCarilli & Walterj 
|2013| >. The Table 2 by |Carilli & Walter| j2013^ gives the average 
CO line ratios for SMGs as Lco{ 4 - 3 )/^co(i-o) = 0.46 and 
^co{ 5 - 4 )/-^co(r-o) = 0.39. Using these ratios, the measured 
CO line luminosities L©o( 4 - 3 ) ±© 0 ( 5 - 4 ) ™ Tablecan be 

translated to ±oo{i-o) of 1 -^ 0.2 x 10 ^^ and 1.1 ± 0.2 x 10 ^^ 

K km s“^ pc“^, respectively. Taking CO(4-3) line luminosity 
(which requires a smaller correction to CO(l-O) line luminosity) 
and a “ULIRG” conversion factor of O. 8 M 0 [Kkms~^ 
we derive a total molecular gas mass of 1.4 ± 0.2 x 10^^M©. Us¬ 
ing a Galactic conversion factor of aco = AT£r2/±oo(i-o) of 
~ 4M© [K kms“^ pc“^]yields a 5 times larger total molecular 
gas mass, 6.8 ± 0.8 x lO^^M©. 

Since we have spectral measurements of multiple CO lines, 
[C II] line, nearly fully mapped SED, and spatially resolved dust 
continuum distribution, we should be able to probe the gas prop¬ 
erties of AzTEC-1 beyond simply adopting a highly uncertain and 
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somewhat arbitrary “average” calibration. In this section, we ex¬ 
plore several different methods for estimating gas mass, including 
dynamical mass analysis, radiative transfer analysis, and dust con¬ 
tinuum measurements, in order to obtain a better handle on the gas 
mass and excitation conditions. 


is highly uncertain because of the unknown and small inclination 
angle, and unfortunately this approach does not offer much useful 
insight into the total gas mass and the CO-to-H 2 conversion factor 
for COSMOS AzTEC-1. 


4.1.1 Gas Mass from Dynamical Mass 


There is a growing awareness that the CO-to-H 2 conversion fac¬ 
tor is not a single value but a quantity dependent on several dif¬ 
ferent factors, such as metallicity, density, temperature, and non- 
gravitational pressure (see a review by |Carilli & Walter||2013| l. 
Since CO is a highly optically thick transition, metallicity is im¬ 
portant mainly for low metallicity systems and should not be an 
important factor for SMGs - they are selected by their bright dust 
continuum and are often bright in CO, [C II], and other transi¬ 
tions that require near solar abundance of metals. The main rea¬ 
son why the “ULIRG” conversion factor is often favored for SMGs 
is that the observed line width may include significant other pres¬ 
sure contributions such as stellar mass and stronger turbulence, as 
well as a higher radiation field resulting from a high density of 
young stars, similar to the central kpc of nearby ULIRGs. |Downes| 
|& Solomon| ( |1998) derived the ULIRG conversion factor, ~5 times 
smaller than the MW value on average, by constructing a dy¬ 
namical model of spatially resolved CO emission in 10 nearby 
ULIRGs and by running a full radiative transfer calculation, also 
taking into account the mass contributions by the new and old stars. 
Because of poor spatial resolution of nearly all existing molec¬ 
ular line imaging data, few examples of such a dynamical mass 
analysis exist for SMGs. In one of the best studies of high red- 
shift SMGs using this approach, |Hodge et al.H20l^ have derived 
OLCO ~ 1-1 ± O. 6 M 0 [Kkms“^ pc“ ]“ for the 2 = 4.05 SMG 
GN20, which has a rather large (14 ± 4 kpc diameter) CO disk 
with a dynamical mass of Mdyn = 5.4 ± 2.4 x stel¬ 

lar mass of M* = 2.3 ± 2.4 x IO^Mq, and a total gas mass of 
Mh2 = 1.8 ±0.7 X 1O“M0. 


We do not have a spatially resolved CO map of COSMOS 
AzTEC-1, but adopting its 890 /rm source size (also the optical 
source siz^ of 0.3” (2.1 kpc) as the diameter of the CO emitting 
molecular gas disk and 1/2 of the observed FWZI [C II] line width 
of 219 km s“^ (see § 4.2 1 as the rotation velocity Vc x {sin i) 
where i is the disk inclination, a dynamical mass can be com¬ 
puted as Mdyn = l.l(sm i)~^ x 1 O^°M 0 . To reconcile this value 
with its stellar mass (1.5 ± 0.2 x 1 O^^M 0 , Smolcic et al.|2011 


or 4.4 X 1O^^M0 according to our GRASIL model) as well as 
the molecular gas mass derived using the ULIRG conversion factor 
(Mh 2 ~ 1 O^^M 0 ), this disk has to be highly inclined, i < 12 °. 
The gas disk traced in [C II] line in the 2 = 4.76 SMG ALESS 73.1 
is 2.2 times larger than the continuum ( |De Breuck et al.|2014| >, and 
the dynamical mass should double if the gas disk emitting [C II] 
in AzTEC-1 is also twice as large. However, the required inclina¬ 
tion changes only slightly, i < 17°. This situation is similar to 
the nearby ULIRG Mrk 231, which has a 450 pc radius molecular 
gas disk with an inclination of i ~ 10° ( [Bryant & Scoville|1996| 
[Downes & Solomon|1998^ . This nearly face-on geometry also of¬ 
fers a natural explanation as why COSMOS AzTEC-1 is excep¬ 
tionally bright in the optical bands. However, its dynamical mass 


4.1.2 Optically Thin CO: a Low Mass Limit 

Approaching the conversion factor problem from an entirely theo¬ 
retical point of view, a lower limit to the total gas mass can be de¬ 
rived by considering an optically thin CO case. [Bryant & Scoville[ 
( [1996[ have offered a rather detailed derivations of gas mass deter¬ 
mination based on spatially resolved CO observations and dynami¬ 
cal modeling in their effort to determine the gas mass for Mrk 231. 
Their derivation shows that the optically thin limit conversion fac¬ 
tor should be aco.min = t).2Q{Xco , where Xco is 
the CO abundance. The minimum gas mass for COSMOS AzTEC- 
1 derived from this relation is ~ 3 x 1 O^°M 0 if the gas excitation 
is typical of SMGs. The minimum gas mass required increases to 
Mh 2 ~ if the excitation temperature is lower as indi¬ 

cated by the measured CO(4-3)/CO(5-4) line ratio (see below). 
These values can be understood better in the context of gas excita¬ 
tion provided the radiative transfer models discussed below. 


4.1.3 Gas Mass from Radiative Transfer Modeling 


The measured ± 00 ( 4 - 3 )/^co( 5 - 4 ) ^ 0.35 

for COSMOS AzTEC-1 is larger than the average SMG ratio of 


0.46/0.39 = 1.18 reported in the recent review by Carilli & Walter 
( [2013[ and is closer to the average Milky Way GMC value of 2.1 - 
see their Table 2. Although omitted by the Carilli & Walter review, 
the scatter in the SMG line ratios found in the literature is rather 
large because many published CO measurements have low S/N 
ratios Qin contrast, the line ratio for AzTEC-1 is quite secure not 
only because each line is detected with S/N > 7 individually, but 
also because the two CO lines are measured simultaneously using 
the RSR, removing any potentially large systematic uncertainties 
arising from utilizing measurements taken at different times using 
different instruments, often on different telescopes. 

We explore the mass and physical conditions of the gas pro¬ 
ducing the observed CO emission in COSMOS AzTEC-1 by ex¬ 
amining a grid of models covering a range of density n, kinetic 
temperature Tmti, and the CO column density Nco that can repro¬ 
duce the measured CO line intensities and line ratios. If the gas was 
optically thin, then only density, temperature and total number of 
CO molecules are needed to model the observed emission. How¬ 
ever, it is more likely that the CO emission is optically thick, and 
this impacts not only the escape of CO photons from the molecu¬ 
lar gas, but also the excitation of CO column trapping. The optical 
depth of the CO emission is determined by the ratio of the CO col¬ 
umn density to the line width AU, and we have run models for 
three cases in which Nco/NV = 2 X 10^'* cm“2 per km s"\ 
6 X 10^® cm“^ per km s“^, and 2 x 10^^ cm“^ per km s“^. In 
the first case the CO emission in the J = 4 — 3 line is optically 
thin (jco !)• In the second case, the optical depth in this line is 
modest (rco ~ 5), and in the final case the optical depth is large 
{tco ~ 20). The excitation by the Cosmic Microwave Background 


^ This stellar source size is estimated from the spatially resolved HST i- 
band image. Toft et al.|(|2014J report an upper limit of r^^NlR < 2.6 kpc 
based on their analysis of the UltraVista near-IR images. 


° For example, in one of the largest recent studies of multiple CO transi¬ 
tions by |Bothwell et al.|j2013^ , 18 out of 32 (56%) CO line measurements 
used for analyzing the SMG line ratios have S/N < 5 - see their Table 5. 
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Table 2. Summary of the RADEX models with acceptable solutions 


n 

(cm-3) 

'^kin ^COl 

(K) 

-0/4cO4-3 

Icoi-o/Icoa-4, 

IcOi-s/Nco 
(Kkms-® cm2) 

(M 0 ) 


Nco/NV 

= 2 X 10®® cm 

1 ^ per km s ^ (optically thin CO cases) 


1 X 10® 

500 

1.28 

2.30 

6.9 X 10"®® 

1.8 X 10®° 

3 X 10® 

200 

0.79 

1.50 

9.6 X 10"®® 

1.3 X 10®° 

1 X 10"® 

100 

0.36 

0.63 

1.4 X 10“®® 

8.9 X 10® 

3 X 10"® 

50 

0.30 

0.54 

1.4 X 10“®® 

8.9 X 10® 

1 X 10® 

33 

0.31 

0.55 

1.2 X 10"®® 

1.0 X 10®° 

3 X 10® 

28 

0.33 

0.58 

1.0 X 10"®® 

1.2 X 10®° 

1 X 10® 

25 

0.37 

0.68 

8.5 X 10“®® 

1.5 X 10®° 


Nco/NV = 6 X 10®® cm-2 

per km s ^ (modest CO optical depth cases) 

1 X 10^ 

120 

2.6 

4.8 

1.5 X 10-®® 

8 X 10®® 

3 X 10® 

80 

1.7 

2.9 

2.4 X 10-®® 

5 X 10®° 

1 X 10® 

35 

1.4 

2.4 

1.8 X 10-®® 

7 X 10®° 

2 X 10® 

25 

1.3 

2.2 

1.1 X 10-®® 

1 X 10®® 

3 X 10® 

20 

1.3 

2.3 

6.3 X 10-®'® 

2 X 10®® 

1 X 10® 

15 

1.3 

2.3 

5.0 X 10-®® 

2.5 X 10®2 

3 X 10® 



No Solution 




Nco /XV ■- 

= 2 X 10®®" cm 

^ per km s ^ (optically thick CO cases) 


1 X 10® 

85 

2.1 

3.6 

7.9 X 10-®'® 

1.6 X 10®® 

3 X 10^ 

35 

1.7 

3.1 

5.2 X 10-®'® 

2.4 X 10®® 

6 X 10® 

25 

1.5 

2.6 

3.4 X 10-®'® 

3.7 X 10®® 

1 X 10® 

20 

1.4 

2.5 

1.9 X 10-®'® 

6.6 X 10®® 

3 X 10® 



No Solution 




(CMB) with Tcmb ~ 14.6 K at 2 = 4.342 is explicitly included 
in these calculations. 


Models producing acceptable solutions are summarized in Ta¬ 
ble]^ The gas density and temperature that satisfy the observed line 
ratio are given along with the ratios of the CO (1-0) line intensity 
relative to the CO (4-3) and CO (5-4) lines and the integrated in¬ 
tensity in the CO (4—3) line divided by the CO column density (in 
units of K km s“^ per cm“^). If we make the assumption that the 
abundance of CO relative to molecular hydrogen is Xco = 10“^, 
then we can determine the amount of mass required to produce the 
observed L'^o in the CO (4-3) line. The last column gives the re¬ 
quired gas mass for that model. 

At a CO column density per unit line width of Nco/XV = 
2 X 10^^ cm“^ per km s“^ where all CO transitions are optically 
thin, an acceptable solution is found over a broad range of tem¬ 
perature. These optically thin cases require the least amount of gas 
mass to reproduce the observations among the models examined, 
but one still needs at least of order of gas to produce the 

observed CO line emission, unless CO is more abundant than it is 
in the Milky Way. A comparison with the optically thin, low limit¬ 
ing mass calculation based on the derivation by [Bryant & Scoville| 
1 1996 I above suggests that the Tco(i-o) estimates derived using 
the average SMG CO line ratios taken from jCarilli & WalterH2013') 
are 2-3 times too large, or the CO emitting gas has to be cold (< 25 
K). 


For the modest to high optical depth cases with Nco /XV = 
6 X 10^® cm“^ per km s“^ and 2 x 10^^ cm“^ per km s“^, ac¬ 
ceptable solutions exist only for a narrower range of density and 
temperature: 10® ^ n < 10® and 15 < Tkin < 120 K. The re¬ 
quired gas mass is at least 5 x 10®® M© in these cases, and more 
likely the range of mass is (1 — 7) x IO^^Mq (a higher mass for 


lower Tkin and higher Nco / NV), comparable or exceeding the 
estimated total stellar mass. If the mean gas density for the star 
forming gas (the “MC” component) is > 10® cm“® as suggested 
by the GRASIL model (see § |3.3.3| >, then there is a very small range 
of parameter space in density and temperature where an accept¬ 
able solution exists, and these solutions favor cold gas temperature 
(< 25 K) with a large total mass {Mh 2 > (2 — 7) x 
An important caveat for this radiative transfer calculation is its as¬ 
sumption of a single component gas. If more than one phase gas 
with vastly different excitation conditions (e.g., “cirrus” and “MC”) 
contribute significantly to the observed CO line intensities and line 
ratios, then a single component analysis such as presented here may 
lead to a misleading result. While a solution satisfying the observed 
^co{ 4 - 3 )/'^co( 5 - 4 ) ‘'3*' found in a fairly broad range 

of excitation conditions and CO column density, an important merit 
of the models summarized Table|^is the clear prediction they make 
on the CO (1-0) line intensity, that the optically thick cases should 
produce five times stronger CO (1-0) line than the optically thin 
cases (for Tkin ^ 200 K). Therefore, future CO (1-0) line mea¬ 
surements should yield an effective discrimination between these 
limiting cases and the gas properties when combined with these 
measurements. 

We briefly explored using [C II] line intensity to gain further 
constraints on the gas excitation and mass, but we found this even 
more problematic. The critical density for excitation of [C II] line 
is much lower (nH 2 ^ 10® cm“®) than those of the CO transitions 
we measured, making this analysis more susceptible to the likely 
presence of multiple gas phases. Furthermore, the [C II] line can 
originate from both neutral and ionized gas (with critical electron 
density of Ue = 10 — 100 cm“®, see Table 2 by 
|2012[ and discussions below), and we do not have much confidence 
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10 Min S. Yun and Other 


in making the assumption that [C II] and CO lines arise from the 
same gas. 


4.1.4 Gas Mass from Dust Continuum 


The molecular gas mass derived from these high J rotational tran¬ 
sitions are likely lower limits since they may be sub-thermally ex¬ 
cited compared with the CO (1-0) transition. One way to check this 
is to compare the total gas mass derived from the Rayleigh-leans 
(Rl) part of the dust spectrum as proposed by |Scoville et al.| ( [2014| l. 
Given the uncertainties in excitation and conversion factor for CO, 
Scoville et al. have argued that dust mass derived from the Rl part 
of the dust spectrum and adopting a gas-to-dust ratio is more robust 
than an estimate based on high J CO line luminosity. The Eq. 12 
by Scoville et al. can be rewritten as 


Mism 


1.2 X 10^0 Tflj 1 S, 
(1-f 2)4-8 ^ To ^ Ujy 


V 

ityiGHz 


r Dl ■ 


A/ 0 (3) 


where Sv is the observed dust continuum in mjy, Dl is luminos¬ 
ity distance in Gpc, and is the RJ correction factor (see their 
Fig. 2). Using this relation, the 345 GHz continuum measurement 
from the SMA (see § |4.2^ can then be translated to a total “ISM 
mass” of Mism = (7.4 ± 1.3) x IO^A/q for Ta = 25 K and 
(3.6 ± 0.7) X 10“A /0 for Ta = 35 K. As noted by Scoville et 
al., estimating dust temperature from the measured dust peak may 
be a mistake since the observed SED is indicative of luminosity- 
weighted (rather than mass-weighted) measure of dust temperature. 
And the resulting smaller gas mass derived for the higher dust tem¬ 
perature is only a lower limit. The error in gas mass depends only 
linearly with any error in dust temperature, and the resulting ~50 
per cent uncertainty due to poorly constrained dust temperature still 
leads to a better gas mass estimate than the estimates from the CO 
luminosity, which is fraught with a wide range of substantial and 
systematic uncertainties. It is notable that the total gas mass de¬ 
rived from the measured dust continuum is among the largest esti¬ 
mates obtained by the different methods, but it agrees well with the 
estimate by the GRASIL model (see § |3.3.3^ and the gas masses re¬ 
quired to produce the observed CO line ratio in the optically thick 
cases (see Table[^. 


Two interesting outcomes from these analysis deserve addi¬ 
tional comments. Firstly, the gas mass analysis using the dynamical 
mass and stellar mass estimates has revealed that the geometry of 
the gas disk has to be nearly face-on, and consequently the poorly 
constrained dynamical mass prevents us from deriving a meaning¬ 
ful gas mass estimate. This calculation also sheds an interesting 
insight in that such a nearly face-on geometry with minimum dust 
obscuration would naturally explain why the stellar component of 
the host galaxy is seen only modestly obscured in the rest frame 
UV light {Av < 3.5), unlike most other SMGs with similarly high 
IR luminosity > Secondly, the observed CO (4-3) to (5- 

4) line ratio is lower than the typical value for SMGs and is closer 
to the Milky Way value. The radiative transfer models for gas den¬ 
sity and temperature characteristic of the MW star forming dense 
cores (n ~ 10"^ cm“® and T ~ 25 K) require a total gas mass of 
(2 — 4) X lO^^Af 0 , which is more in line with the gas mass esti¬ 
mate from the dust continuum and the gas mass estimate based on 
the MW value for aco (~ 5 x lO^^A/ 0 ). The range of excitation 
conditions that can reproduce the observed line ratios and intensi¬ 
ties are uncomfortably narrow, however, and this may indicate that 
the observed CO lines include significant contributions from more 
than one component of molecular ISM present in this galaxy (e.g., 
[Harris et al.|2010^ , as also suggested by the GRASIL analysis. 

Given the range of gas mass estimates and current star forma¬ 
tion rate, the gas depletion time for COSMOS AzTEC-1 is about 
200 Myr, with about a factor of 2 overall uncertainty. This means 
COSMOS AzTEC-1 will exhaust its gas reserve and will shutoff its 
star formation activity by 2 « 4 even without any negative feed¬ 
back, unless gas continues to flow in at a rate matching the star 
formation rate, M « IO^Mq yr~^. Some gas recycling can extend 
this time by about 50 per cent, but the stellar feedback process is 
expected to do more than compensating for this effect. The stel¬ 
lar mass doubling time r* = A/* /SFR for COSMOS AzTEC-1 is 
also about 200 Myr, and its substantial stellar mass could have been 
plausibly built up entirely during the current episode of starburst 
(see jYun et al.|2012l ). In such a scenario, the starburst activity would 
have started around 2 « 5, ending with a M* > 6 x stel¬ 

lar galaxy with ~ 2 kpc diameter by 2 « 4, similar to the massive 
quiescent galaxies reported byj Whitaker et al.H2013| l, |Straatman'^ 
[^j2014| l, and others. 


4.1.5 Summary of Gas Mass Estimation and Broader 
Implications 


4.2 [C IIJ/FIR Ratio and High Radiation Field 


A detailed review of the process of converting the measured CO 
(4—3) and (5-4) line luminosity to a total gas mass demonstrates 
several assumptions one has to make, particularly when using em¬ 
pirical calibrations. The measured line ratio between these two tran¬ 
sitions for COSMOS AzTEC-1 is larger than the average ratio re¬ 
ported for a sample of SMGs by jCarilli & Walt^j2013| >, and the 
resulting uncertainty in estimating CO (1-0) line luminosity is large 
(a factor of 2 to 3). The CO-to-H 2 conversion factor is also uncer¬ 
tain by at least a factor of 2 or more, but the nominal total molecular 
gas mass based on the CO (4-3) line luminosity is 1.4 x 10^^ A/ 0 . 
The non-LTE radiative transfer calculations have yielded a range 
of acceptable solutions summarized in Table with likely molec¬ 
ular gas mass in the range of (1 — 7) x MY^Mq for the modest to 
high optical depth cases and with a minimum (optically thin) limit 
of ~ 10^°A/ 0 . The total ISM mass estimate based on dust contin¬ 
uum jScoville et al.|20i4l l favors the upper end of these estimates, 
(4 — 7) X while the empirical calibration that yields a to¬ 

tal gas mass of ~ 2 x HY^Mq from the measured CO (4-3) line 
luminosity is at the low end of these different estimates. 


The 158 /im [C II] line is an important coolant of the neutral ISM 
and thus is a bright tracer of star formation in galaxies, typically ac¬ 
counting for 0.1-1 per cent of IR luminosity ( [Madden et al.|19^ 
[Malhotra et al.|2dOT] [Stacey et al.|2010[ l. Because [C II] emission 
can be produced by different gas phases with a wide range of phys¬ 
ical conditions, interpreting [C II] emission is difficult (see a recent 
review by [Goldsmith et al.|2012[ (. A broad correlation is seen be¬ 
tween observed [C II] emission and other tracers of star formation 
(e.g., [Boselli et al.|[200^ |de Looze et al.[[2011[ also see Fig. |^. 
However, interpreting observed [C II] line luminosity in terms of a 
particular physical process, such as a tracer of SFR, is problematic 
because [C II] emission arises from a variety of different excitation 
mechanisms, in both ionized and neutral phase. Also, star forming 
galaxies observed in [C II] show a factor of 100 or more spread 
in the Lycii]/L fir ratio, which is correlated with IR luminosity 
and dust temperature (so-called “[C II] deficiency”, [Malhotra et al.[ 
[19971 ILuhman et al.|1998> . 

The measured [CII] line luminosity of COSMOS AzTEC-1 is 
significantly higher than those of the Great Observatories All-sky 
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L'FIR [Lq] 


Figure 7. as a function of far-IR luminosity Lpjfi. Empty cir¬ 

cles are the GOALS sample of local LIRGs and ULIRGs (|Diaz-Santos et| 
|al.|2013) . Empty squares and crosses are [C II] measurements for z > 1 


sources from the literature jWalter et al.||2009| |Hailey-Dunsheath et al. 

2010 

IStacey et al.|2010| |Wagg et al.|2010||Cox et al.|2011| |Valtchanov 

et al. 

2011||Gallerani et 

al. 120121 ISwinbank et al.|20121|Venemans et al. 

2012 

[Walter et al.|2012 

ICamiani et al.|20131 (George et al.|2013|lWang 

et al. 

20131 IMagdis et a 

.|2014| [Rawle et al.|2014|[Riechers et al.|2013[ 

Willott et al.||2013| |Riechers et al.||2014| |Brisbin et al.||2015| |Gullberg 

et al. 

|2015| (Schaerer et al. 20l5k - crosses are strongly lensed sources 


while squai'es may also be lensed sources. COSMOS AzTEC-1 is shown 
as a large filled circle and extends the trend of “[C II] deficiency” to 
Lfir ^ IO^^Lq. Filled squares are high redshift sources with spatially 


resolved [C II] and FIR distribution: COSMOS AzTEC-3 jRiechers et al. 

2014(, HDF 850.1 [Neri et al. 

2014k BR1202-0725 A&B ICarniani et al. 

2013), ALESS 73.ipe Breuc 

c et al.|2014), and four z = Q QSOs imaged 


using ALMA jWang et al.|2013^ T A typical error bar is shown on the bottom 
right comer. 


LIRG Survey (GOALS) sample of 241 luminous infrared galax¬ 
ies studied by |Diaz-Santos et al.| ( |2013l l using the Herschel Space 
Observatory (see Figure |7J. Also shown are a collection of [C II] 
line sources at z > 1 from the literature, and they extends the ob¬ 
served broad correlation to Lfir > 10^^Lq. The [C II]/FIR ratio, 
shown in Figurej^ reveals that the measured L[c ii] /T fir ratio of 
6.5 X 10“^ for AzTEC-1 is among the lowest measured and extends 
the [C II] deficiency to Lfir IO^^Lq. 

While most of the GOALS sample LIRGs and ULIRGs form 
a broad trend with a decreasing L[c ii\ / Lfir ratio with increasing 
Lfir, the z > 1[C II] sources show a much larger scatter. These 
high redshift systems simply being a scaled up versions of the lo¬ 
cal star forming galaxies is a commonly offered explanation (e.g., 
[Stacey et al.|20ld{|Brisbin et aLpOlS) . At least 1/2 of the sources 
detected in [C II] thus far are strongly lensed systems (shown as 
crosses in Fig.ISl found by the South Pole Telescope (SPT; [^ 


berg et al.[2015|( and Herschel iCox et al. 

2011| [Valtchanov et 

al. 

20111 Riechers et al.l|2013| (Magdis et a 

.|[2014| Rawle et al. 

2014 1 , and many should fall along the local LIRG/ULIRG relation 


when corrected for magnification, as shown by a detailed study of 
a z = 2.013 lensed [C II] source by [Schaerer et al.| ( |2015^ . De¬ 
termining whether the larger scatter associated with the remaining 



Figure 8 . [C II]/FIR ratio as a function of the FIR luminosity. Many of 
the 2 > 1 sources are lensed, and their luminosities are not corrected for 
lensing because the magnification factor is not always known. All symbols 
are identical to those in Figure|7] A typical error bar is shown on the bottom 
left comer. 


2 > 1 sources can be accounted by lensing will require future de¬ 
tailed follow-up studies of the individual sources. 

Along with the 2 = 5.3 SMG COSMOS AzTEC-3 ([Riech- 


ers et al.[[M4l ) and the 2 = 5.2 SMG HDE 850.1 ( [Neri et al. 


2014i, AzTEC-1 has the one of the smallest L[c ii]/L fir ra¬ 
tio in Figure clustered together with 7 IR luminous galaxies 
hosting an optical QSO at 2 > 4 imaged in [C II] and contin¬ 
uum by ALMA (shown as filled squares). One possible explana¬ 
tion for their extremely low L[c ii]/L fir ratio is the elevated 
dust-obscured AGN contribution to the FIR luminosity, but this 
explanation is not supported by systematic studies of large sam¬ 
ples of IR luminous galaxies conducted using Herschel. By exam¬ 
ining the GOALS LIRGs with and without AGN activity (identified 
through Spitzer mid-IR spectroscopy), [DIaz-Santos et aq j2013[ l 
have shown that the correlation between L\c ii\ and Lfir and the 
[C II] deficiency is an intrinsic property of the star formation, and 
there is no need to invoke AGN activity to explain it (at least at lev¬ 
els of > 10“^). A study of 154 intermediate redshift 

(( 2 ) ~ 0.15) 24 /r-m-selected galaxies by |Magdis et al. pQ13^ and 
a study of 130 mid- and far-IR selected galaxies by|Sargsyan et al.| 
( [2014^ also drew a similar conclusion. 

The compact source sizes of AzTEC-1 and other sources re¬ 
vealed by high resolution continuum imaging using the SMA and 
ALMA suggests the high intensity of the infrared radiation field 
may offer an important clue to the [C II] deficiency. Possible ex¬ 
planations for the [C II] deficiency include: (1) self-absorption; (2) 
saturation of the [C II] line due to high gas column density; (3) 
decreased photoelectric heating in high UV radiation field; and (4) 
high dust-to-gas opacity caused by an increase of the average ion¬ 
ization parameter (see reviews by|Malhotra et al.|2001||DIaz-Santos| 
[et al.[2013| (. Citing a clear trend for LIRGs with deeper 9.7 p,m sili¬ 
cate strengths, higher mid-IR luminosity surface densities (Emir), 
smaller fractions of extended emission, and higher specific star for¬ 
mation rates (SSFRs) to display a greater [C II] deficiency, | Dfaz-[ 
[Santos et al.[ ^2013^ have concluded that the dust responsible for 
these correlations must be directly linked to the process driving the 
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[Vkpc=] 


Figure 9. [C IIJ/FIR ratio as a function of FIR surface density. Only the 
sources with a resolved mid-IR or far-IR sizes are included. All symbols 
are identical to those in Figure|7] A typical error bar is shown on the bottom 
left corner. 


observed [C II] deficiency. They have also found the correlation 
becoming much tighter when the FIR luminosity is normalized by 
mid-IR source size (i.e., surface density Ef/a), independent of the 
nature of the powering source. As shown in Figure COSMOS 
AzTEC-1 and other high redshift [C II] sources with spatially re¬ 
solved continuum sizes follow the same tight correlation defined by 
the local LIRGs and extend this correlation by one order of mag¬ 
nitude larger in Ef/h- In addition to extending this correlation to 
a higher luminosity density, this comparison also further supports 
the proposed scenario that the compactness of the active region 
and the resulting higher intensity of the infrared radiation field dic¬ 
tates the [C II] deficiency. In an earlier modeling study using the 
spectral synthesis code CLOUDY, [Gracia-Carpio et al.| ( |201 f| have 
shown that all far-IR fine structure lines, regardless of their origin 
in the ionized or neutral phase of the ISM, show a deficit with in¬ 
creasing Lfir/Mh2 ratio, and they further conclude that this de¬ 
ficiency is driven by the increased ionization parameter. This is an 
extremely interesting prediction that should be tested further using 
future observations of other far-IR fine structure lines. Determin¬ 
ing through a high resolution imaging study whether the unlensed 
Lfir > [C II] sources with Lyc ii]/L rir ^ 10“® fol¬ 

low the narrow trend seen in Figurej^is another important test for 
this [C II] deficiency scenario. 


5 CONCLUSIONS 

We report the first successful spectroscopic redshift determination 
of COSMOS AzTEC-1 obtained with a clear detection of the red¬ 
shift CO (4-3) and CO (5-4) lines using the Redshift Search Re¬ 
ceiver on the Large Millimeter Telescope and the confirmation of 
the CO redshift through the detection of the redshifted 158 pm 
[C II] line using the Submillimeter Array. Utilizing the newly mea¬ 
sured redshift and CO and [C II] line intensities, we have explored 
the gas mass and physical conditions of the gas fueling the enor¬ 
mous luminosity associated with this a = 4.342 SMG. 

The RSR spectrum of COSMOS AzTEC-1 (Figure[T]( has two 


emission lines clearly above the noise level (^ la) at 86.31 GHz 
and 107.87 GHz, and they are identified as redshifted CO (4-3) 
and (5-4) lines at z = 4.3420 ± 0.0004, respectively. This con¬ 
clusion is supported by its over-all SED as well as by the radio- 
millimetric spectral index analysis by |Carilli & Yun| ( |1999| l. A de¬ 
tailed discussion of the unique redshift determination is presented 
in the Appendix section. The derived CO redshift of j* = 4.3420 
is slightly lower than the photometric redshift derived by|Smolcic| 
| et al-lpOl l| l using the rest frame UV and optical photometry data 
and is outside the failed previous blind CO searches by|Smolcic et| 
[aT HIOTT} and |Iono et al.| ( [20T2l l. This successful redshift determi¬ 
nation after nearly 10 years of effort demonstrates the power of the 
ultra-wideband spectroscopic capability of the RSR on the Large 
Millimeter Telescope. The redshifted 492 GHz [C I] line is not de¬ 
tected (S[c i]/Sco(i- 3 ) ^ 0.45), but this upper limit is still in 
line with the measured [C I] line strengths in other high redshift 
galaxies i Walter et al.|201 1 


[C I]/ScO(3-2) 


■ 0.3). 


The new CO redshift for COSMOS AzTEC-1 is verified by 
the detection of redshifted [C II] line at 335.8 GHz using the SMA. 
The bright [C II] line is detected with S/N ~ 15, and a higher 
spectral resolution clearly shows that the line is asymmetric. The 
cause of this asymmetry is not known yet, but this explains the 
slightly lower redshift determined for the [C II] line. Although the 
derived [C II] line luminosity of Lyc ii] = 7.8 x 10® Lq is re¬ 
markably high, it is only 0.04 per cent of the total IR luminosity, 
making COSMOS AzTEC-1 one of the most [CII] deficient objects 
known. We show that AzTEC COSMOS-1 and other high redshift 
[C II] sources with a spatially resolved source size extend the tight 
trend seen between the [C II]/FIR ratio as a function of FIR surface 
density among the IR-bright galaxies by |DIaz-Santos et al.H2013} 
by more than an order of magnitude. This result lends further sup¬ 
port for the explanation that the higher intensity of the IR radiation 
field and the resulting increased ionization parameter are likely re¬ 
sponsible for the “[C II] deficiency” seen among luminous infrared 
starburst galaxies. 

Our modeling of the observed spectral energy distribution us¬ 


ing a modified black body model, starburst SED models by Ef- 


stathiou et al.| ( |2000| l, and the GRASIL SED code ( [Silva et ah 


199811 produces a consistent estimate of the IR luminosity (Lir 
(1.4-1.6)xlO^®Lo andiF/H = (0.9-1.2) x lO^^L©). The esti¬ 
mated star formation rate from the IR luminosity (1600—1700 M© 
yr“^) is slightly larger than the model-based SFR (880 and 1320 
M© yr“^ for the Efstathiou and GRASIL model, respectively). The 
model SFR and total stellar mass estimates depend on the adopted 
star formation history which is intrinsically more uncertain. The 
best fit GRASIL model constrained by the observed luminosity and 
the shape of the UV-to-radio SED further suggests of an intense, 
compact starburst (r* « 0.1 kpc) heavily obscured {Ay > 200 for 
the molecular clouds and Ay = 3.5 for the “cirrus” component) 
by a massive, compact gas cloud (Mgas = 3.6 ± 0.6 x lO^^M©, 
rgas ~ 1 kpc). 

The total molecular gas mass was derived from the mea¬ 
sured CO (4-3) and CO (5-4) lines and 345 GHz continuum us¬ 
ing several different methods, specifically addressing the uncer¬ 
tainties associated with each method. A grid search for non-LTE 
radiative transfer models that match the observed CO line inten¬ 
sity and line ratio yields acceptable solutions over a wide range of 
gas temperature and density with a minimum (optically thin) limit 
of ~ 10 ^® M©. However, plausible models (modest to high opti¬ 
cal depth) require a narrow range of gas temperature (T ~ 20- 
35 K) for densities n > 10^“® cm“®, requiring gas masses of 
Mh 2 = (1 — 7) X lO^^M©, Conventional methods of comput- 
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ing molecular gas mass from the observed CO line intensities are 
subject to very large uncertainties in translating these high J transi¬ 
tions to the intensity of the CO (1-0) line, as well as to the similarly 
uncertain “aco” conversion factor. The empirical calibration that 
yields a total gas mass of ~ 2 x from the measured CO 

(4-3) line luminosity is on the low end of these different estimates. 
The total ISM mass derived from the 345 GHz continuum ( |Scoville| 
|et al.|2014|l is near the top of the mass range derived by the other 
methods: Mjsm = (7.4 ± 1.3) x IO^Mq for Td = 25 K and 
(3.6 ±0.7) X for Td = 35 K. Future measurements of the 

CO (1-0) transition should remove the uncertainty associated with 
the translation of the higher J lines and offer a useful constraint on 
the CO optical depth. 

Our dynamical mass analysis shows that the gas disk in COS¬ 
MOS AzTEC-1 has to be nearly face-on in order for the derived 
dynamical mass to be consistent with the minimum possible com¬ 
bined gas and stellar masses. This offers a natural explanation for 
the bright, compact stellar light distribution visible in the rest frame 
UV band HST images, similar to the situation in the local ULIRG 
Mrk 231. The same analysis also suggests extremely high opacity 
{Av > 200) for most other viewing angles, as seen in many other 
high redshift SMGs. 

Among the 15 brightest AzTEC sources identified by the 
AzTEC/JCMT survey of the COSMOS field and located with a bet¬ 
ter than 1" positional accuracy using the SMA observations, COS¬ 
MOS AzTEC-1 is only the second object with a secure spectro¬ 
scopic redshift, after the z = 5.3 COSMOS AzTEC-3 jRiechers~et| 
|al.|2010^ . Advent of the RSR on LMT and other similar broadband 
spectrometer systems on modern telescopes with a large collect¬ 
ing area (e.g., ALMA) is finally making accurate determination of 
redshifts of these distant, optically faint galaxies possible. In addi¬ 
tion to yielding redshifts, these CO spectroscopic surveys can also 
yield information on total gas masses and dynamical masses along 
with the excitation conditions of the gas fueling the rapid growths 
of these young, massive galaxies. A complete RSR survey of these 
COSMOS AzTEC sources has started at the LMT, and we should 
soon be able to gain an unbiased view of the redshift distribution 
and total molecular gas mass contents of these and other SMGs. 
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APPENDIX A: REDSHIFT DETERMINATION FROM AN 
RSR SPECTRUM 

Al Template Cross-correlation Analysis 

The simultaneous frequency coverage of the RSR between 73 and 
111 GHz means at least one CO transition falls within the spectral 
coverage at all redshifts except for a narrow redshift range between 
0.58 < 2 < 1.08, and two or more CO or [C I] transitions fall 
within the RSR spectral range at z ^ 3.15 (see Fig. |Al| l. A vari¬ 
ety of fainter molecular transitions from less abundant species such 
as HCN, HCO+, HNC, CS, CN, HC 3 N, and H 2 O have also been 
detected in nearby and distant galaxies (see a review by|Carilli &| 
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Figure A2. Template cross-correlation amplitude of COSMOS AzTEC-1 
RSR spectrum in S/N unit. A zoomed in details of the most significant peak 
(S/N = 9.0) at z = 4.342 is shown in the inset, and it clearly shows that 
the CO lines are clearly resolved spectrally. 



Figure Al. Observed frequencies of redshifted CO and [C I] line transitions 
falling within the RSR frequency coverage range (73 to 111 GHz). At least 
one CO line should appear in the RSR spectrum at all redshifts except for a 
nan'ow redshift range of 0.58 < 2 : < 1.08. Two or more CO & [C I] lines 
should appear simultaneously within the RSR spectrum at 2 > 3.15. 


Walter||2Q13| and references therein). As first introduced by |Yun| 


et al. 1 2QQ7| |, a cross-correlation analysis is a powerful method to 


derive the redshift information from such a broadband spectrum, 
even when many of the lines are not individually detected with a 
good S/N ratio. A cross-correlation product (^{z) can be derived 
as a function of redshift 2 : from the observed spectrum S{iy) and 
the model spectral template M{p^z) as 


C(^) = J S{p)M{p,z)W{u)du. 

The Doppler shifted model spectral template M(t/, z) is derived as 

r(i^+Ay/2){l+z) 

M{p,z)— / T{v)dp' 

J (f^ —A^'/2)(1 + 2) 

where T{v') is the rest frame template spectrum and At/ is the 
RSR channel width. The weight function W{v) represents the rel¬ 
ative strength of different molecular transitions, and an empirical 
composite spectrum based on observed relative line strengths for 
high redshift sources (e.g., |Spilker et al.|2014| l is adopted for the 
analysis presented here. 

The number of spectral lines contributing to the model spectral 
template M{u, z) increases with redshift as the total frequency cov¬ 
erage of the RSR in the rest frame grows as 38(1-1-2:) GHz. As a re¬ 
sult, the noise in the cross-correlation amplitude ^( 2 ) increases ac¬ 
cordingly with redshift, and interpreting the raw cross-correlation 
amplitude is not straightforward. Also, since many of the molecular 
transitions occurring in the millimetre and sub-millimetre bands are 
rotational transitions with only slightly different rotation constants, 
the distribution of line transitions is highly clumped in the spectral 
domain, further complicating the situation. Therefore, rather than 
interpreting the raw cross-correlation amplitude ^( 2 ) for the red¬ 
shift analysis, we compute a ‘"S/N ratio” of ((( 2 ) for a quantitative 
analysis of acceptable redshift solutions. The “noise” in each red¬ 
shift bin is estimated by randomly shuffling the input RSR spectmm 


10,000 times, and the derived cross-correlation amplitude ((( 2 ) is 
converted to a histogram of S/N ratio as shown in Figure [A2] 


A2 Determination of A Unique Redshift Solution 

The histogram of the template cross-correlation amplitude for the 
COSMOS AzTEC-1 RSR spectrum in Figure [A2| has the highest 
peak with S/N = 9 at 2 = 4.342, but other peaks with an appar¬ 
ent S/N > 5 are also seen. Since the cross-correlation analysis is 
sensitive to all real signal, the two spectral line features detected 
near 86 GHz and 108 GHz in Figure [T] eac/r produce a series of 
S/N=5-l peaks that corresponds to different rotational transitions 
of CO at 2 ~ 0, 1, 2, & 3, in addition to the strongest peak re¬ 
sulting from the two CO lines at 2 = 4.342. The presence of two 
distinct lines for AzTEC-1 in the RSR spectrum (Fig.f^ rules out 
all single line identifications at 2 < 3.15 (see Eig. |AT^ , and the 
2 = 4.324 solution remains as the only plausible interpretation 
with S/N > 5. A small but non-zero possibility of two unrelated 
CO sources at 2 < 3.15 along the same line of sight still remains 
(e.g., [Zavala et al.|2015^ , but we can rule out this scenario by using 
other tests. 

Photometric redshift constraints can be extremely helpful for 
determining the likely redshift identification ( |Yun et al.|2007| l. Ex¬ 
ploiting the well known radio-IR correlation among star forming 
galaxies and the strong positive and negative fc-corrections at the 
radio and millimetre wavelengths, the radio-millimetric spectral in¬ 
dex technique ( jCarilli & Yun|[T99^ in particular is a simple but 
remarkably powerful method that requires just two broadband pho¬ 
tometry measurements. As shown in Eigure [A3] the product of the 
probability distribution for radio-millimetric photometric redshift 
p(z) and 5'Ai?[^(2)] (shown in Fig. |A2[ l effectively removes all 
2 < 3 scenarios and nicely isolates the 2 = 4.342 solution. The ef¬ 
fectiveness of the radio-millimetric photometric redshift method is 
further demonstrated by the fact that a powerful constraint against 
low redshift solutions can be derived even when only a good ra¬ 
dio upper limit is used - all low redshift solutions can be re¬ 
jected equally well by treating the ~ 4 (t radio photometry point 
of AzTEC-1 only as an upper limit. 
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Figure A3. A plot of the template cross-coiTelation amplitude of COSMOS 
AzTEC-1 RSR spectrum in S/N unit with the redshift constraint from the 
radio-millimetric spectral index technique jCarilli & Yun|1999) utilizing 
just the AzTEC 1.1mm and the VLA 1.4 GHz photometry. 

























